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We revisit in detail tlie non-mean-field ground-state phase diagram for a binary mixture of spin- 
1 Bose-Einstein condensates including quantum fluctuations. The non-commuting terms in the 
spin-dependent Hamiltonian under single spatial mode approximation make it difficult to obtain 
exact eigenstates. Utilizing the spin z-component conservation and the total spin angular momen- 
tum conservation, we numerically derive the information of the building blocks and evaluate von 
Neumann entropy to quantify the ground states. The mean-field phase boundaries are found to 
remain largely intact, yet the ground states show fragmented and entangled behaviors within large 
parameter spaces of interspecies spin-exchange and singlet-pairing interactions. 
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I. INTRODUCTION 

Ultracold atomic quantum gases with spin degrees of 
freedom provide exeeptionally clean and idealized test- 
ing beds for studying quantum magnetism [l|, Q . Optical 
trapping from ac Stark shifts of ofF-resonant laser fields 
are capable of equal confinement for all atomic pseudo- 
spin components which facilitates research into exciting 
spinor physics with atomic quantum gases. As is often 
employed in studying a trapped Bose gas, when treat- 
ing the condensate, we first take a mean-field (MF) ap- 
proximation assuming only one eigenvalue of the single- 
particle density matrix is macroscopic, being of order N . 
Second quantization is then limited to the condensate 
mode. Such a simple scenario already allows for many 
interesting quantum many body phenomena^]. 

Two popular atomic species often employed in ex- 
perimental research on spinor Bose-Einstein condensates 
(BECs) are ^^Rb and ^SNa atoms. Within each 

species, their interactions are dominated by the density- 
dependent interaction in comparison to the much weaker 
spin-dependent interactions. As a result, single-spatial 
mode approximation (SMA), whereby the spatial depen- 
dence of the condensate wave function is determined in- 
dependent of the spin degrees of freedom, was introduced 
and remains reasonable as long as the number of atoms 
is not too large [l^. Within the MF approximation, the 
ground state of a spinor Bose-Einstein condensate (BEG) 
is found to be ferromagnetic, polar, or cyclic phases, etc, 
determined by the spin-dependent interactions and the 
total (hyperfine-)spin F of the atom. Further theoret- 
ical work armed with full quantum calculations reveal 
interesting many-body states d, 0, H, [HI, beyond the 
scope of those from MF approximations. For spin-1 con- 
densates, the exact eigenstates will contain paired spin 
singlets 0, 0, [HI , which become more complex for higher 
spin condensates. For example, the spin-2 case involves 
spin singlets which can be formed by either two or three 
atoms 0, Hi- A general procedure exists for more de- 
tailed information of the building blocks of eigenstates 
determined by their associated generating functions [l^ . 



More generally, we can always resort to the means of 
numerics to diagonalize the ground-state single-particle 
density matrix, which then reveals fragmented ground 
states if more than one eigenvalues are being of order 

Several groups have recently studied spinor conden- 
sate mixtures |17H24j . which also display non-MF fea- 
tures, such as anomalous quantum fluctuations for each 
spin components and quantum entangled ground states. 
This is first discovered in spinor condensates with more 
than one orbitals, for instance, the case of pseudo spin- 
1/2 condensates |14| - [l6| for which Kuklov and Svistunov 
14| predicted that in the ground states all atoms will 
have to condense into two orthogonal spatial orbitals or 
more due to the conservation of the total spin. This 
could result in a condensate ground state being a maxi- 
mally entangled many-body state. Shi et al. replaced the 
two orbitals with two different atomic species, a ground 
state with entangled order parameter follows [16| . Under 
MF approximation, we have previously elaborated the 
ground-state phase diagram for a condensate mixture of 
two spin-1 BECs [l^. The interesting phases are named 
appropriately as FF, AA, PP, CC, and MM phases, dis- 
tinguishing different structures and interaction parame- 
ter spaces. Furthermore, we provide many beyond MF 
results based on a full quantum spin-dependent Hamil- 
tonian [l9| . which contains non-commuting terms, for- 
bidding a simple derivation of the exact eigenstates. For 
two special cases, commutations are restored among the 
generally non-commuting terms. First, when the inter- 
species singlet pairing interaction is ignored (7 = 0), the 
Hamiltonian is then simply composed of three operators 
which obey the angular momentum algebra. Making use 
of the eigenstates of single spin-1 condensates, we directly 
construct the eigenstates of a binary spin-1 mixture using 
the angular momentum coupling representation. Second, 
when the interspecies anti-ferromagnetic spin-exchange 
interaction is strong enough, the ground state will be 
forced to develop entanglement between the two species 
[igj , a result consistent with what is discovered in a spin-1 
condensate placed inside a double well [1^. Other inter- 
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esting features are discussed for 7 = revealing fragmen- 
tation and quantum entanglement [13, Hi], [2^ . 

II. THE MODEL HAMILTONIAN 

In this revisit we hope to understand quantum entan- 
glement between two spin-1 condensates when both in- 
terspecies spin-exchange and singlet-pairing interactions 
are present. Our study is based on the same model sys- 
tem of a binary mixture of spin-1 condensates confined 
in optical traps. The corresponding field operators that 
annihilates a boson of species 1 and species 2 at posi- 
tion r are described respectively by ^MF(r) and <&Mj7(r), 
where Mp = —1,0,1 denoting the three Zeeman hyper- 
fine states. The SMA is adopted for each of the two 
species, employing two spatial mode functions ip{v) and 
(/)(r) respectively, and the field operators are expanded 
as *Mi.(i") = ctMpijir) and l>Afj.(r) = foAfF0(r), with 
and respectively the annihilation operators for 
an atom in the spin component Mp. In the absence of 
external magnetic field, the spin-dependent Hamiltonian 
for a binary mixture of spin-1 condensates then becomes 
the following 

Hs = iCi/3i (L?-2iVi) +ic2/?2(^2-2iV2) 

+ iCi2/3il-i2 + ^Ci270l2012, (1) 
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under the SMA [l^. The interaction coefficients are 
Ci = /dr|^(r)|4, C2 = /dr|0(r)|4, and C12 = 
/ c?r|i/)(r)p|(/)(r)p. (3i {[i^) is intra-species spin-exchange 
interaction parameter of species 1 (2). /3 and 7 denote 
inter-species spin-exchange and singlet-pairing interac- 
tion parameters, respectively. The singlet pairing op- 
erator becomes G)|2 = o\V_i — a\b\ -I- aLi^i, and two 
angular momentum like operators Li = "Y^ij '^t^ij^j ^i^*^ 
L2 = Sii_^FijSj obey the usual angular momentum al- 
gebra [6l.l26j. They commute with atom number opera- 
tors Ni = ajdi and N2 = X^i ^l^i- In the above, 
denotes the component of the spin-1 matrix F. 

III. GROUND-STATE PHASE DIAGRAM 

As presented in the earlier study of [l^, for the spin- 
dependent Hamiltonian of Eq. ([T]), the first three terms 
commute with each other, but they do not commute with 
the fourth term. This shows the ground state determined 
will depend on the interaction parameters. As a result, 
we resorted to the special cases of no interspecies singlet- 
pairing interaction (7 = 0) and Ci/3i = C2P2 = C12P/2 
hjl. For the first case of 7 = 0, it has already at- 
tracted much attention due to the appearance of frag- 
mented ground states and the associated entanglement 
between two species and exotic atomic number fluctua- 
tions [ii-ii|. 



In this study, we will discuss the general case of the 
full spin-dependent Hamiltonian of Eq. ([T]). Whenever 
the ground state depends on the interaction parameters, 
we have to perform a full quantum calculation numeri- 
cally, usually this amounts to a full exact numerical di- 
agonalization for each atom numbers iVi and N2. Be- 
fore discussing the numerical results, we want to point 
out that there still exist two conserved quantities: the 
total spin angular momentum and its ^-component, as 
= [Ll + L2)^ commutes with the spin-dependent 
Hamiltonian. As a result, we can elaborate spin struc- 
tures from building blocks derived by generating function 
of the maximum spin states \l,lz = I), where we have 
used quantum numbers I and 1^ to denote the common 
eigenstates of the angular momentum operators and 
L,. 

We recall the suitable generating function Gg{x,y,z) 
for a binary mixture of two spin-1 condensates derived 
earlier in Ref. [T9| . From this generating function, we 
have figured out all six building blocks for constructing 
the eigcnstate \l,l), which is given by 

where the six building blocks are 
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and Ui ,Vi, and Wi satisfy the constrains 

Ui + 2U2 + Wi + W2 ~ Ni, 
Vi + 2V2 + Wi + W2 = N2, (4) 
Ui + Vi + W2 ~ I, 

and additionally W2 ~ 0,1. Spin states \l,lz 7^ I) of 
other magnetization can be constructed by simply ap- 
plying {LS'"' on the state \l,l) as L''~'''\1,1) oc \l,lz 7^ I) 
(un-normalized) . 

In most numerical studies, we assume each species con- 
tain 100 atoms (Ni = N2 = N = 100). Due to the SO (3) 
symmetry of our model, we restrict the Hilbert space into 
the subspace of zero magnetization with Iz ^ [131 • 
Fig. [U we summarize the extensive numerical results. 
Since the total spin angular momentum is conserved, we 
can use the eigenvalue of the operator to distinguish 
different phase, which is then accompanied by the infor- 
mation of the building blocks. A total of three constrains 
exist for the allowed values of u^, Vi, and Wi (i = 1,2). 
Only three are needed for a solution, we choose the three 
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as U2 , V2 and wi , which are determined numericahy from 
evaluating the associated expectation values of Ag Aq , 
B^^''^B^^\ and f [i^'^^'^f respectively. 

From extensive numerical results we construct the 
ground-state phase diagram as shown in Fig. [T] Per- 
haps not surprisingly, it is almost the same as the MF 
approximation studied in Ref. [l^. Each phase is then 
labeled the same as before [l^, albeit that the meanings 
can be different due to the non-commuting operators in 
the spin-dependent Hamiltonian Eq. ([1]). In Tab. HI we 
summarize the properties for the four special phases: FF, 
AA, PP, and CC. The remaining MM phase still denotes 
the phase whose parameters evolve continuously across 
the phase boundaries. 













FF 


= 2N{2N+1) 


= 


= 
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AA 
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~ N{N + 2) 


PP 


= 


~ A'(iV-|- 1) 


~ iV(iV-f 1) 


~ 


CCl 


= N{N + 1) 


~ N{N + 1) 


~ 


~ A(iV + 2) 


CC2 


= 


~ N{N + 1) 


~ A(iV + 1) 


~ N{N + 2) 


CCS 


= A(iV + l) 


~ 


~ A(iV + 1) 


~ A(A + 2) 



TABLE I: The expectation values for the special operators in 
the ground state within different phases. 



IV. ENTANGLED GROUND STATES 

To quantify entanglement between the two species, we 
numerically computed the von Neumann entropy S{pi) = 
— Tr(/5i log2jv+i Pi), where pi = Tr2P is the reduced den- 
sity matrix resulting from partial tracing of the ground- 
state density matrix p over the basis of species 2. The 
amount of entanglement is then shown as density plots 
over the phase diagram of Fig. [1] with a legend shown 
on the right, the black (white) color refer to low (high) 
entanglement. 

In the absence of interspecies singlet-pairing inter- 
action (7 = 0), the spin-dependent Hamiltonian of 
Eq. ([T]) contains only three operators commuting with 
each other. As a result we can use quantum numbers ^i, 
I2, I, and Iz to quantify its eigenstate \h,l2il,lz), with 
(L?) = kih + 1), (Ll) = hih + 1), {L'^) - l{l + 1) and 
(Lz) — Iz- In the ground states, we have I = li + I2 
for ferromagnetic interspecies spin-exchange interaction 
(P < 0); and I = |^i — Z2I for anti-ferromagnetic inter- 
species spin-exchange interaction (/3 > 0), while the value 
of h and I2 are determined by the three interaction pa- 
rameters. 

In Fig. [21 we display the expectation values of the 
intra- and interspecies single-pairing number operators 
in the ground state, divided by their corresponding max- 
imum values shown in Tab. HI In addition, we show the 
total spin angular momentum and the von Neumann en- 
tropy. First, when < 0, /32 < 0), irrespective of the 
interspecies spin-exchange interaction, atoms in the same 
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FIG. 1: (Color online). The ground-state phase diagram and 
the corresponding von Neumann entropy distribution at fixed 
values of Ci/3i and C2/32. Blue solid lines denote continuous 
phase transition boundaries. Red solid lines denote discon- 
tinuous phase transition boundaries between two phases with 
fully determined total spin angular momentum I. The green 
solid line denotes the discontinuous phase transition boundary 
between the phase CCS and the MM phase. The black dash- 
dotted lines correspond to C127 ~ 0, which serve as guides for 
the eye. The von Neumann entropy of the ground states are 
painted by gray scale density plots, where black (white) color 
refer to low (high) entanglement respectively. The three sub- 
plots denote fixed intra-specie spin exchange interaction pa- 
rameters of (Ci^i,C2/32)/|Ci^i| — (a) (-1,-2) ; (b) (1,2); 
and (c) ( — 1,2). The red dashed lines are four lines con- 
necting the points Oaa(— S2, 6), Occ(32, 6), C'ff(32, —6) and 
Opp(— 32, —6) in the parameter space of (C127, Ci2/3)/|Ci/3i|. 




species will not pair into a singlet, but atoms in different 
species will pair into singlets with (Fg^'^'^-'^rQ^'"'^-') close to 
reaching its maximum values N{N + 2) when the inter- 
species spin-exchange interaction is anti-ferromagnetic. 
The corresponding total spin angular momentum is 
equal to its maximum value 2N{2N + 1) with a rela- 
tively low von Neumann entropy when the interspecies 
spin-exchange interaction is ferromagnetic. The total 
spin angular momentum is equal to with the von 
Neumann entropy S{pi) = 1 for anti- ferromagnetic inter- 
species spin-exchange interaction. Second, when (/3i > 
0, (32 > 0), in the two limits of large ferromagnetic or anti- 
ferromagnetic interspecies spin-exchange interaction, the 
ground state is the same as above in the previous case. 
In the other limit with low interspecies spin-exchange in- 
teraction, atoms in the same species tend to pair into 
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singlets giving rise to no entanglement between the two 
species. This implies the ground state can be written 

as a product state: Z-^/'^{A^o'>'')'^/^{B^^'>^)'^/^\vSic). For 
the remaining phase, we call it the MM phase, which can 
show higher (lower) entanglement compared to the FF 
phase between the two species when /? > (/3 < 0). For 
the final case when < 0,(32 > 0), the ground state 
show similar properties as that in the second case. 
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FIG. 2: (Color online). LEFT: The ground-state normalized 
expectation values of intra- and interspecies singlet-pairing 
number operators A'^^^'^ A'^^^ /N{N + 1), B^^^"^ B^^^ /N{N + 1), 
p(i,i)tp(i,i)/^(^ + 2), denoted by blue solid, red dash-dot, 
and black dash-dot-dot lines, respectively. RIGHT: The nor- 
malized total spin angular momentum /2N{2N + 1) and 
von Neumann entropy of the ground state, denoted by blue 
solid and red dash-dot lines, respectively. The three subplots 
denote zero interspecies singlet-pairing interaction (7 = 0) 
and fixed intra-species spin-exchange interaction parameters 
(Ci/3i,C2/32)/|Ci/?i| =: (a) (-f,-2) ; (b) (1,2); and (c) 
(-1,2). 

The most attractive phase when 7 = is the entangled 
ground state denoted by ipj^j^^ p^ . 

= J2 (-)^-™|iV,m)i®|iV,-m)2,(5) 

which show high entanglement with S{pi) ~ 1 be- 
tween the two species. As demonstrated in Fig. [T] 
by numerical calculations, however, there remain other 
phases which show larger entanglement between the two 
species. This shows that the state tpj^j^ is not a max- 
imal entangled state, in contrast to previously studied 
case of two pseudo spin- 1/2 condensates [l^. This is 
not a surprise [2^. Due to the redundant degrees of 



freedom in the spin-1 case, the total spin angular mo- 
mentum of each species can take other values besides 
the largest value of A^. To demonstrate the entangle- 
ment between the two species, we show their correspond- 
ing expectation values and von Neumann entropy along 
four lines connecting the points Oaa(— 32, 6), Occ(32, 6), 
Off(32,— 6) and Opp(— 32,— 6) in the parameter space 
of (C127, Ci2/3)/|Ci/3i|. The four lines are marked as red 
dashed lines in the Fig. [TJa). 

In Fig. [21 we illustrate the ground-state proper- 
ties of two ferromagnetic condensates with intraspecics 
spin-exchange interactions at (Ci/3i, C2;82)/|Ci/3i| = 
(— 1, —2). First of all, we consider the AA phase. When 
7 = 0, the spin-dependent Hamiltonian contains three 
operators commuting with each other, and the ground 
state can be expressed as Vaa ^O'" large enough anti- 
ferromagnetic interspecies spin-exchange interaction, and 
show high entanglement between the two species. When 
7 7^ 0, although the fourth term of the spin-dependent 
Hamiltonian does not commute with the other three, we 
find that the ground state not only show similar expecta- 
tion values of the operators, it also contain similar entan- 
glement between the two species, over a large area in the 
phase diagram demonstrated in Fig. [IJa). For 7 < 0, 
irrespective of its value, the ground state falls into the 
AA phase. For 7 > 0, the ground state is still classi- 
fied as the AA phase, as long as C127 does not exceed 
a critical value, which increases in proportion to the in- 
terspecies spin-exchange interaction parameter Ci2/3. In 
the first column of Fig. [3l we evaluate the properties of 
the AA phase, where the ground-state expectation value 
of intra- and interspecies singlet-pairing number opera- 
tors are close to 0, 0, and N{N + 2), respectively. The 
total spin angular momentum is exactly equal to 0, and 
the von Neumann entropy is close to 1. 

When interspecies singlet-pairing interaction exceeds 
a critical value, the ground state changes into the MM 
phase. As long as C127 > 0, it tries to decrease the in- 
terspecies singlet-pairing interaction. In the first column 
of Fig. [31 we follow the line of OaaOcc and illustrate 
the phase transition from the AA phase into the MM 
phase. In the MM phase, as long as we increase C127 
accordingly, atoms in different species will continuously 
avoid to pair into singlets, while atoms in species 1 will 
try to pair into singlets. Atoms in species 2 first try to 
pair into singlets and then avoid to pair. Meanwhile, the 
total spin angular momentum first increases and then de- 
creases. For a relatively large area of the MM phase, we 
find the two species show high entanglement compared 
to the AA phase. 

As interspecies singlet-pairing interaction is increased, 
the ground state will fall into the CCl phase, where its 
total spin angular momentum will be equal to N(N + 
1), and atoms in the species 1 (2) will pair (not pair) 
into singlets. At the same time the expectation value of 
To^'^'^rl,^'^' will be near to its maximum N{N + 2). The 
ground-state von Neumann entropy in the CCl phase 
remains at a low value. 
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FIG. 3: (Color online). TOP: The ground-state normalized expectation values of intra- and interspecies singlet-pairing number 
operators A^^^ A')^'> /N(N + 1), B^q ^'^B^^^/N{N + 1), f '^'^'"''f ^^'^'/^(^^ + 2), denoted by blue solid, red dash-dot, and black 
dash-dot-dot lines, respectively. BOTTOM: The normalized total spin angular momentum /2N{2N -f 1) and Von Neumann 
entropy of the ground state, denoted by blue solid and red dash-dot lines, respectively. From left to right, we illustrate 
the corresponding expectation values and Von Neumann entropy along four direct lines connecting four points Oaa(— 32,6), 
Occ(32, 6), Off(32, — 6) and Opp(— 32, — 6) in the parameter space of (Ci2/3, Ci27)/|Ci/3i |. The intraspecies spin-exchange 
interactions are fixed at (Ci/3i, C2/32)/|Ci/3i| = (—1, —2). The four lines are marked as red dashed lines in the Fig. [TJa). 
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FIG. 4: (Color online). The same as that in Fig. [3] but with different intraspecies spin-exchange interactions fixed at 
(Ci/?i,C2/32)/|Ci/3i| = (l,2). 



Going along the line OccOff, with decreasing in- 
terspecies spin-exchange interaction, the ground state 
changes from the CCl phase, to the MM phase, and fi- 
nally to the FF phase. Atoms of species 1 will become 
unpaired continuously, while the total spin angular mo- 
mentum increases from N{N -|- 1) in the CCl phase to 
its maximum 2N(2N + 1) in the FF phase. The ground- 
state entanglement between the two species in the CCl 
and the FF phase are almost at the same level. While in 
the MM phase, the entropy first increases to near 1 and 
then decreases. 

We then follow the line OffOpp. With decreasing 
interspecies singlet-pairing interaction, the ground state 
covers the FF, MM, and PP phases successively. In the 
PP phase, atoms in the same/different species will all try 



to pair into singlets. The ground-state expectations val- 
ues for intra- and interspecies singlet-pairing number op- 
erators will reach close to their corresponding maximum, 
meanwhile the two species show higher entanglement. In 
the MM phase, the expectations values for operators or 
the von Neumann entropy change continuously to con- 
nect the FF phase and the PP phase. 

Lastly, we consider the line OppOaa- As long as 
the interspecies spin-exchange interaction increases, the 
ground state will be changed continuously from the PP 
to the MM, and to the A A phase. From the fourth col- 
umn of Fig. [3l we find that in the whole line of OppOaa, 
the two species show a relatively high entanglement with 
the von Neumann entropy S{pi) remaining higher than 
1. Especially so in the MM phase, where the highest 
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FIG. 5: (Color online). The same as that in Fig. O but with different intraspecies spin-exchange interactions, which are fixed 
at (Ci/3i,C2/32)/|Ci/3i| = (-1,2). 



entropy reaches 1.4568, which is close to the maximum 
entropy of log^jy+^HN + 1){N + 2)/2) ~ 1.6116. 

In Figs. 13] and [51 we illustrate the two 
other cases with intraspecies spin-exchange interac- 
tions (Ci/3i, C2/32)/|C'i/3i| fixed respectively at (1,2) and 
(—1,2). We find that the ground-state show similar prop- 
erties to that in the two ferromagnetic condensates shown 
in Fig. |3l The only difference is for the CC2 or CCS 
phase. In the CC2 phase, the ground state expecta- 
tion values for both intra- and interspecies singlet-pairing 
number operators are close to their corresponding maxi- 
mum N{N+1), N{N+1), and N{N + 2). While in the 
CCS phase, they arc close to 0, iV(iV-l-l), and iV(iV-|-2), 
respectively. The total spin angular momentum of the 
ground state is equal to in the CC2 phase, and N{N+1) 
in the CCS phase. Meanwhile, in the CC2 phase the von 
Neumann entropy is less than 1 but larger than that in 
the FF phase, while in the CCS phase, it is close to the 
value in the FF phase. 

Before conclusion, we hope to stress that the maximal 
entangled state in this system is given by 

V'ME = ^-'/'(f(''''^)'^|vac). (6) 

It is the eigenstate or the ground state (if 7 < 0) of the 
7-term in the Hamiltonian of the Eq. ([1}, which means 
the maximal entangled state V'me is the eigenstate of 
two spin-1 condensates with only interspecies spin-singlet 
pairing interaction (/3i = ^2 ~ P ~ and 7 7^ 0), with 
the corresponding eigenvalue Ci2jN{N + 2)/6 [29| . 



more carefully, going beyond the MF approximation. 
When there exists no interspecies singlet-pairing inter- 
action, the spin-dependent Hamiltonian contains three 
operators commuting with each other. In this special 
case, the most interesting phase is the AA phase, where 
two species show high entanglement. When interspecies 
singlet-pairing interaction is turned on, the added opera- 
tors do not commute with the previous three ones, which 
forbids us from obtaining exact eigenstates analytically 
for the model spin system. In this study, we perform 
full quantum diagonalization to find the ground states 
numerically. To quantify the ground states, we work 
out the building blocks to construct the maximum spin 
states, which can be used rightfully to discuss entangle- 
ment scales between the two species. We have evaluated 
the associated ground-state von Neumann entropy. After 
detail calculations, we find that the AA phase can per- 
sist for large areas of the parameter space for interspecies 
spin-exchange and singlet-pairing interactions. In addi- 
tion, there is another interesting phase: the PP phase, 
which show similar level of entanglement between the 
two species. What's more, we find that AA phase is not 
the maximum entangled state. The ground state with 
highest entanglement we found lies in the MM phase. 
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